Naval  Research  Laboratory 

Washington,  DC  20375-5320 


NRL/MR/5580-18-9767 


Node  Ranking  Tool  -  NoRT 


Ira  S.  Moskowitz 
Jonathan  W.  Decker 

Information  Management  &  Decision  Architectures  Branch 
Information  Technology  Division 


Rebekah  H.  Kang 

University  of  Virginia 
Charlottesville,  VA 


March  23,  2018 


DISTRIBUTION  STATEMENT  A:  Approved  for  public  release.  Distribution  is  unlimited. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including 
suggestions  for  reducing  this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway, 
Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of 
information  if  it  does  not  display  a  currently  valid  0MB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY) 
23-03-2018 


2.  REPORT  TYPE 

Memorandum  Report 


3.  DATES  COVERED  (From  -  To) 
October  2015  -  September  2017 


4.  TITLE  AND  SUBTITLE 


Node  Ranking  Tool  -  NoRT 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 

62235N 


6.  AUTHOR(S) 

Ira  S.  Moskowitz,  Jonathan  W.  Decker  and  Rebekah  H.  Kang* 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 

IT-235-018 


5f.  WORK  UNIT  NUMBER 

6804 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

U.S.  Naval  Research  Laboratory 
4555  Overlook  Avenue,  SW 
Washington,  DC  20375-5320 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 


NRL/MR/5580-18-0039 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Office  of  Naval  Research 

One  Liberty  Center 

875  North  Randolph  Street,  Suite  1425 

Arlington,  VA  22203-1995 


10.  SPONSOR  /  MONITOR’S  ACRONYM(S) 

ONR 


1 1 .  SPONSOR  /  MONITOR’S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

DISTRIBUTION  STATEMENT  A:  Approved  for  public  release.  Distribution  is  unlimited. 


13.  SUPPLEMENTARY  NOTES 

*  University  of  Virginia,  Charlottesville,  VA 


author’s  email:  ira.moskowitz@nrl.navy.mil 


14.  ABSTRACT 

This  paper  gives  a  description  of  the  Node  Ranking  Tool  (NoRT)  that  was  developed  as  part  of  the  Applied  Network  Science  6.2  base  program 
work  unit  at  NRL,  Code  5580.  We  explain  the  theory  of  NoRT  and  how  to  use  it. 


15.  SUBJECT  TERMS 

TOPSIS,  Social  Network,  Sensor  Network,  Centrality,  Diffusion,  Disease,  Virus,  Expectation,  Pandemic,  Closeness,  Graph,  Degree,  Spectrum 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT 

18.  NUMBER 

OF  PAGES 

19a.  NAME  OF  RESPONSIBLE  PERSON 

Ira  S.  Moskowitz 

a.  REPORT 

Unclassified 

Unlimited 

b.  ABSTRACT 

Unclassified 

Unlimited 

c.  THIS  PAGE 

Unclassified 

Unlimited 

SAR 

23 

19b.  TELEPHONE  NUMBER  (include  area 
code) 

(202)  404-7930 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39.18 


This  page  intentionally  left  blank. 


11 


1 


U.S.  NAVAL 
LresearckJ 

LABORATORY 


Node  Ranking  Tool — NoRT 

Ira  S.  Moskowitz*,  Jonathan  Decker^,  Rebekah  Kang^ 

*  Information  Management  &  Decision  Architectures  Branch 
Code  5580,  Naval  Research  Laboratory 
Washington,  DC  20375 

§  Information  Management  &  Decision  Architectures  Branch 
Code  5581,  Naval  Research  Laboratory 
Washington,  DC  20375 

^  NRL  Code  5580,  NREIP  Summer  Intern  2017 
University  of  Virginia 
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Abstract 

This  paper  describes  the  Node  Ranking  Tool  (NoRT)  that  was  developed  as  part  of  the  Applied  Network  Science 
6.2  base  program  work  unit  at  NRL.  We  explain  the  theory  of  NoRT  and  how  to  use  it. 

Index  Terms 

TOPSIS,  Social  Network,  Sensor  Network,  Centrality,  Diffusion,  Disease,  Virus,  Expectation,  Pandemic,  Close¬ 
ness,  Graph,  Degree,  Spectrum. 


I.  Introduction 


THis  paper  gives  a  description  of  the  Node  Ranking  Tool  (NoRT)  that  was  developed  as  part  of  the  Applied 
Network  Science  6.2  base  program  work  unit  at  NRL,  Code  5580.  We  explain  the  theory  of  NoRT  and  how 
to  use  it.  Examples  are  supplied.  We  first  explain  centrality,  then  we  explore  some  applications,  and  then  we  turn 
to  a  discussion  of  NoRT  itself. 

We  wish  to  determine  which  nodes  in  a  network  are  “strong”  and  which  are  “weak”.  We  model  a  network  as 
a  connected  undirected  graph  G  with  a  set  of  nodes  N  =  {oi,  ...,0^}.  G  has  as  an  adjacency  matrix  A,  which  is 
a  symmetric  matrix  that  has  i,j  entry  1,  if  and  only  if  there  is  a  link  between  Ui  and  oj,  and  is  zero  otherwise. 
We  denote  the  link,  if  it  exists,  between  Ui  and  oj  as  lij  =  1.  If  there  is  no  link,  then  lij  =  0.  The  set  of  {kj}  is 
denoted  by  L.  Note  that  all  lij  have  a  constant  weight  of  1  or  0,  and  there  are  no  self-loops.  Hence,  as  noted  A 
is  symmetric  with  zeros  down  the  diagonal  and  O’s  or  I’s  as  the  off  diagonal  entries,  where  Aij  =  weight(ti_j). 

In  this  paper,  the  network  and  the  graph  are  taken  as  identical  entities.  We  envision  an  attack  on  the  network 
in  the  form  of  nodes  becoming  “compromised”.  Once  a  node  is  compromised,  it  acts  as  an  “agent”  ll25l  of  the 
attack  and  attempts  to  compromise  other  nodes.  We  wish  to  rank  the  network  nodes  in  terms  of  their  ability  to 
compromise  other  parts  of  the  network,  or  to  be  compromised  themselves. 

In  our  view,  a  network  being  compromised  via  the  spread  of  a  “disease”  @|.  Once  a  disease  is  in  a  node(s),  we 
are  concerned  with  how  the  disease  spreads  to  the  rest  of  the  network.  Of  interest  is  the  asymptotic  behavior  of  the 
disease,  particularly  in  the  way  it  affects  the  other  nodes  in  the  long  term.  That  is,  do  nodes  become  healthy  once 
they  are  infected?  Do  all  the  nodes  eventually  become  infected  resulting  in  a  pandemic?  Do  certain  parts  of  the 
network  become  infected,  while  others  are  left  untouched?  We  may  also  be  concerned  with  how  a  disease  spreads 
over  a  short  period  of  time. 
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We  use  a  combination  of  various  forms  of  graph  centrality  to  determine  the  likelihood  of  a  node  to  compromise 
other  parts  of  the  graph,  or  to  become  compromised  itself.  Furthermore,  we  also  consider  the  asset  criticality 
ifTTl  of  the  graph  nodes.  That  is,  some  nodes  are  more  important  than  others  in  terms  of  becoming  infected,  or 
spreading  disease.  Our  study  of  disease  propagation  must  therefore  be  fine-tuned  to  consider  asset  criticality,  which 
is  important  if  we  wish  to  plan  defensive  or  offensive  measures  involving  certain  specific  nodes. 

We  balance  all  of  the  above  concepts  by  using  the  ideas  as  given  in  171  [121  HSl  (which  discusses  Technique  for 
Order  of  Preference  by  Similarity  to  Ideal  Solution  (TOPSIS)  which  is  a  Multi- Attribute  Decision  Making  (MADM) 
technique),  node  importance,  and  asset  criticality.  These  concepts  will  be  made  clear  in  the  rest  of  the  document. 

1-A.  What  is  NoRT? 

NoRT  is  an  application  of  the  TOPSIS  method  where  there  are  four  criteria,  consisting  of  the  different  node 
centralities  discussed  below.  Furthermore,  NoRT  has  a  fifth  optional  criteria  called  asset  criticality  which  allows  us 
to  individually  denote  nodes  that  are  more  important  than  others. 

NoRT  is  a  self-contained  SageMath  script  run  from  the  command  line.  The  code  is  included  in  this  paper. 

Experiments  have  been  run  to  validate  the  theoretical  basis  of  using  NoRT.  Extensions  to  the  present  NoRT  code 
are  easily  done. 


II.  Node  Centrality 

There  are  many  different  measure  of  node  centrality  and  importance  ll26l  Ch.  7].  We  concentrate  on  four  in 
NoRT.  They  are  (1)  Degree  Centrality  (DC),  (2)  Closeness  Centrality  (CC),  (3)  Betweenness  Centrality  (BC),  and 
(4)  Eigenvector  Centrality  (EC). 

Every  connected  graph  has  a  natural  distance  function  for  its  nodes.  We  say  that  the  distance  between  Vi  and  Vj 
is  given  by 

d{ui,Uj)  :=  shortest  (in  terms  of  the  number  of  links)  path  length  between  Oi  and  Uj.  (1) 

Given  and  uj,  we  say  that  a  path  between  Oi  and  Uj  is  a  geodesic  iff  the  length  of  that  path  is 
Note  that  a  geodesic  between  Ui  and  Uj  need  not  be  unique.  We  denote  the  set  of  geodesics  between  Oi  and  i/j 
as  G{vi,Uj).  We  see  that  for  the  3E  graph  in  Eig.  1  that  G(l,4)  =  {(1-0,  0-4),  (1-3,  3-4)},  a  set  with  two  paths. 


Eig.  1:  3E  Graph 


However,  G(0,  3)  is  made  up  of  three  geodesics  paths  {(0-1, 1-3),  (0-2,  2-3),  (0-4, 4-3)}. 
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11- A.  Degree  Centrality  (DC) 

DC  is  simply  the  normalizecQ  degree  (the  number  of  links  into  a  node)  of  a  node.  For  node  i>i  we  have 

DC{vi)  :=  ^  ^i,j  ■  (2) 

j 

Since  our  graph  is  connected,  DC  can  never  he  zero,  and  the  most  links  a  node  can  have  is  m  —  1.  Therefore 
DC{vi)  ranges  from  to  1,  where  m  is  the  cardinality  of  the  set  of  nodes. 


Fig.  2:  Star  graph  7  nodes 

In  Fig.  I^we  see  that  every  node  is  connected  to  vi,  which  has  maximal  DC  =  1,  and  all  the  other  nodes  have 
DC  =  1/6. 

DC  is  a  reliable  measure  for  how  quickly  one  node  can  send  “activity”  out  over  a  link  and  compromise  other 
nodes,  especially  if  the  compromising  ability  weakens  over  time,  or  is  probabilisitic  in  nature.  That  is,  one  can 
view  DC  as  measuring  how  a  diffusive  wave  spreads  over  time  from  a  node,  with  its  amplitude  diminishing  over 
time/hops. 


11-B.  Closeness  Centrality  (CC) 

For  node  Oi,  we  consider  d{ai,i'j),  take  its  inverse,  and  normalize  by  m  —  1.  Again,  we  note  that  the 
normalization  is  optional;  however,  we  use  the  same  algorithms  as  SageMath  |Q.  In  general,  our  centrality  measures 
may  differ  from  others  in  the  literature  by  constants.  This  does  not  affect  node  rankings  so  is  insignificant  to  us. 
The  closeness  centrality  (CC)  of  node  Ui  is  given  as 


CC{u,) 


m  —  1 


'when  we  introduce  the  TOPSIS  we  will  see  that  how  attribute  values  are  normalized  is  irrelevant. 


(3) 
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Thus  for  Fig.  we  have  that  CC{vi)  =  1  and  CC{vj)  =  2-b+i-i  ~  3  ~  ^■ 

Centrality  closeness  basically  tells  us  that  a  node  that  is  closer  to  the  other  nodes  has  a  larger  CC  than  one  that 
is  further  away.  A  good  discussion  of  CC  can  he  found  in  |[2^  7.6].  Closeness  centrality  gives  us  a  measure  of 
how  quickly  a  disease  may  spread  from  one  node  to  the  rest  of  the  nodes. 


11-C.  Betweenness  Centrality  (BC) 

Another  important  measure  of  centrality  is  BC.  Given  node  t'i,  we  consider  two  subsets  of  all  the  geodesics  of 
G.  The  first  is  simply  G{vj,Vk),  as  defined  above,  and  the  other  is  G\i,^{i'j,Vk),  which  is  the  subset  of  G{i'j,i'k) 
consisting  of  all  geodesics  between  Uj  and  Vk,  provided  that  Vi  is  an  interior  (not  end)  point  of  the  geodesic.  We 
define  fhe  Befweenness  cenfralify  as: 


BC{n,) 


fm-l\  ^  ^  \\Gl,{nj,Uk)\\ 

V  2  )  ||G(.,,.,)|| 


(4) 


A  high  BC  value  is  indicative  of  being  a  very  critical  node  in  terms  of  being  an  intermediary  for  information 
transmission. 


Looking  at  Fig.  we  see  that  node  1  is  interior  to  every  geodesic  between  the  other  nodes.  We  only  have  one 
geodesic  between  any  two  nodes  2,...,7,  since  a  geodesic  does  not  care  about  direction,  this  gives  us  15  geodesics 
without  node  1  as  an  endpoint.  However,  each  one  of  these  geodesics  goes  through  node  1.  Hence 

BG{1)  =  Q  •  15  =  1 

which  makes  the  strange  normalization  term  clear.  But  on  the  other  hand  BG{j)  =  0,  j  =  2, ...,  7. 

However,  node  j,  j  =  2, ...,  7  is  never  interior  to  any  geodesic  in  G,  so  G\y\vi,  v^)  =  %\i,k  arbitrary,  j  =  2, ...,  7. 
Therefore,  BG{j)  =  0,j  =  2,  ...7  in  Fig.  2. 


Now  looking  at  Fig.  we  calculate  BG{1).  The  only  subset  of  the  geodesics  that  contains  node  1  as  an  interior 
point  is  G(0,  3).  In  this  subset  there  is  only  one  geodesic  that  has  node  1  as  an  interior  point.  Since  (^2^)  =1/6, 

we  have  that  BG{1)  =  1/18.  Appealing  to  symmetry,  we  also  see  that  BG{2)  =  BG{4)  =  1.  Now  consider 
BG{0),  the  subsets  that  have  node  0  as  interior  points  are  G(l,4),  G(l,2)  and  G(2,4).  Since  each  of  these 
subsets  has  two  elements,  with  only  one  element  per  subset  having  node  0  as  an  interior  point,  we  have  that 
BG{0)  =  1/6  •  (1/2  +  1/2  +  1/2)  =  1/4.  Again,  by  symmetry  BG{3)  =  1/4  also. 


II-D.  Eigenvector  Centrality  (EC) 

Our  final,  and  extremely  important,  measure  of  centrality  is  EC.  Note  there  are  many  other  measures  of  centrality. 
We  feel  that  the  four  in  this  document  suffice  for  mosf  nefwork  vulnerabilify  analyses.  Additional  measures  can 
easily  be  added  info  NoRT,  if  need  be. 

The  adjacency  matrix  A  of  G  is  obviously  non-negative.  Therefore  we  have  the  following  famous  result  lIT^  p. 
536,  Thm,l]  and  ||T9l  p.  543,  Thm,l]  . 

Theorem  2.1:  (Perron-Frobenius)  A  has  a  real  eigenvalue  Xmax  greater  than  or  equal  to  the  magnitude  of  any  of 
its  other  eigenvalue^  There  is  a  nonnegativ^  eigenvector  Ua  associated  with  Xmax- 

To  settle  the  choice  of  normalization  once  and  for  all,  we  use  the  Gould  vector  lllQl,g,  as  the  unique  representative 
of  Ua  with  Euclidean  length  1.  We  use  the  notation  Ui  to  denote  the  unit  vector  that  is  0  in  all  components,  except 
for  the  ith  component  which  is  1.  We  define  Il2l.ll26l  Ch.  7]  the  Eigenvalue  centrality  (EC)  of  node  Vi,  to  simply 
be  the  i  component  of  94, 


^The  spectral  radius  of  A  is  also  called  Xmax- 
^Of  course,  Ua  is  unique  only  up  to  a  positive  scalar. 


EC{vi)  :=  Q-Ui  . 


(5) 
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The  spectral  radius  is  used  as  a  phase  transition  point  for  various  models  in  virology.  The  relative  size  of  the 
components  of  the  Gould  vector  influences  the  spectral  radius.  What  this  implies  is  that  the  most  effective  place 
to  put  a  disease,  in  terms  of  infiltration,  is  in  the  node  that  has  the  largest  component  value  of  the  Gould  index, 
see  |[28l  |2^  and  Shield  value  Q.  EC  is  also  a  key  part  of  many  page  rank  algorithms. 

For  Fig.  12  we  have  that  Amax  =  VG,  and  9  =  ^  '  I/a/S,  I/a/S,  I/a/G,  1/v^;  1/VG}-  Not  surprisingly 

the  node  (rjin  the  middle  has  the  highest  EC  value. 


II-E.  Interesting  Example  1:  the  Kite 


Consider  Fig.  Which  is  the  correct  metric  of  centrality?  Fet  us  make  a  table  of  the  different  centrality  measures 
with  the  node  ranking  next  to  the  values. 

We  see  that  the  node  ranking/importance  is  very  much  dependent  upon  which  centrality  measure  we  use. 
Therefore,  following  Q  and  Q,  we  use  TOPSIS  to  evaluate  our  decisions  about  the  relative  importance  of  the 
nodes  in  a  graph.  Of  course,  in  this  example,  DC  and  EC  display  the  same  ranking,  but  for  other  graphs  that  need 
not  be  the  case.  However,  the  raw  values  of  these  ranks  affects  the  TOPSIS  process.  For  example,  in  Fig.  4  the  six 
unique  (highest  to  lowest)  EC  rankings  are  nodes  1,5,6,4,3,2,  and  the  three  DC  rankings  are  1,5,(3,4,6),2.  Where 
0  means  they  are  the  same  ranking. 
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Node  in  Rank  Order 

DC*9  Value 

7 

6 

4 

5 

5 

5 

1 

4 

2 

4 

3 

3 

6 

3 

8 

3 

9 

2 

10 

1 

Node  in  Rank  Order 

CC  Value 

4 

0.64 

5 

0.64 

7 

0.6 

8 

0.6 

1 

0.53 

2 

0.53 

3 

0.5 

6 

0.5 

9 

0.43 

10 

0.31 

Node  in  Rank  Order 

BC  Value 

8 

0.39 

4 

0.23 

5 

0.23 

9 

0.22 

7 

0.1 

1 

0.02 

2 

0.02 

3 

0 

6 

0 

10 

0 

Node  in  Rank  Order 

EC  Value 

7 

0.48 

4 

0.4 

5 

0.4 

1 

0.35 

2 

0.35 

3 

0.29 

6 

0.29 

8 

0.2 

9 

0.05 

10 

0.01 

TABLE  I: 

(f-V) 


Rankings  against  Four  Different  Centrality  Measures  for  the  Kite  Fig.  3.  Note  that  BC  is  normalized  hy 
instead  of  to  agree  with  [[hi  (this  has  no  effect  on  the  rankings). 


II-F.  Questions 

If  we  were  to  “attack”  the  network,  which  node  would  we  attack  first?  If  we  were  to  defend  the  network,  where 
would  we  put  the  most  defensive  firewalls?  If  we  were  to  start  removing  links  to  thwart  an  attack,  at  which  node(s) 
would  we  remove  the  links?  If  some  nodes  are  more  critical  than  others,  how  would  that  influence  the  above?  We 
will  show  later  in  this  paper  how  TOPSIS  can  answer  these  questions. 

III.  Social  Networks 

We  quote  from  Il24l 

Social  networks  are  a  core  part  of  modern  information  sharing  and  communication  and  distinctly 
related  to  Internet  media  consumption.  As  a  result,  understanding  information  flow  within  these  networks 
is  of  great  interest  to  both  industry  and  government  EOl.  Social  networks  are  used  by  roughly  20%  of 
the  people  in  the  United  States  to  make  decisions  ll22l.  which  is  a  direct  indication  of  the  importance 
of  social  networks  such  as  Facebook,  Twitter,  etc.  Within  this  context,  the  transmission  of  information 
across  the  social  network  topology  can  provide  indications  of  events  and  situational  awareness  of  all 
sorts,  e.g.  terrorist  activities,  marketing  penetration,  social  consciousness,  etc.  ||T1.  Thus  researchers  have 
applied  many  techniques  to  understand  information  transmission/diffusion  in  the  characterization  of  this 
phenomenon.  The  diffusion  of  innovation  over  a  network  is  one  of  the  original  reasons  for  studying 
networks.  Furthermore,  the  spread  of  disease  among  a  population  has  been  studied  for  centuries  ifTTl. 

The  basis  for  many  of  the  social  network  analysis  tools  that  seek  to  understand  or  predict  diffusion  is 
the  concept  of  centrality  ll2^  [lH  @1 .  In  most  applications  of  social  network  analysis,  centrality  typically 
is  used  to  identify  a  node/vertex  that  is  the  most  influential  or  having  the  most  importance.  In  the 
case  of  flow,  spread,  or  diffusion,  this  can  be  taken  as  a  super  spreader  |l2l.  For  brevity,  we  begin  by 
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Fig.  4:  EC  and  DC  rank  differ 


highlighting  the  importance  and  utility  of  centrality  to  social  network  analysis.  A  thorough  discussion  of 
centrality  measures,  as  well  as  some  of  their  limitations,  can  he  found  in  Freeman’s  work  JH.  Freeman 
correctly  pointed  out  that  centrality  must  he  considered  in  the  context  of  the  network  structure/topology 
and  offers  centrality  as  control,  centrality  as  independence,  or  centrality  as  activity.  Relative  to  social 
network  information  diffusion,  and  our  work  here,  we  consider  centrality  in  all  three  contexts  and  focus 
on  a  significant  consideration  in  its  use  within  that  domain. 

Well-studied  metrics  of  centrality  ll2^  Ch.  7]  do  not  take  prohahility  into  account.  In  ll24l  some  surprising 
examples  were  shown  emphasizing  that  a  study  of  graph  behavior  based  solely  on  topological  and  degree  properties 
is  incomplete  when  it  comes  to  modeling  virus  or  information  propagation.  We  emphasize  the  spectral  radius  of 
a  connected  graph  because  when  there  is  both  an  infection/information  transmission  and  cure  rate,  the  inverse 
spectral  radius  gives  one  a  phase  transition.  Furthermore,  the  eigenvector  corresponding  to  the  spectral  radius  gives 
us  a  very  important  measure  of  node  centrality.  However,  we  de-emphasize  the  cure  rate,  which  is  consistent  with 
information  transmission  such  as  might  occur  within  a  social  network  (this  paragraph  paraphrased  from  CH). 

In  Il20l  a  detailed  study  of  “influencers”  is  given.  They  also  discuss  the  interesting  inverse  problem  of  determining 
where  a  social  influence  started  and  potentially  controlling  mis-information.  They  discuss  three  information  diffusion 
models  in  social  networks  and  how  the  centrality  measures  of  degree,  betweenness,  closeness,  and  eigenvalue  type 
centrality  must  be  analyzed  in  different  venues  of  social  network  analysis.  In  |[^  Sec.  2.2.2]  they  discuss  how 
different  measure  of  centrality  for  identifying  influencer  nodes  are  discussed  in  the  literature. 

In  |[2T1  different  measures  of  centrality  are  analyzed  with  respect  to  media  coverage  of  actors  in  social  media. 
In  lITSl  Sec.  2.2.4]  measures  of  centrality  similar  to  what  we  use  are  discussed  in  detail  and  compared. 

IV.  Sensor  Networks 

In  ifTH  the  authors  discuss  the  need  for  more  work  on  centrality  measures  for  sensor  networks.  The  same  authors 
show  further  results  with  respect  to  eigenvector  centrality  and  randomly  deployed  sensor  networks  in  lIT^. 

We  note  that  temporal  changes  in  sensor  networks  can  be  studied  using  the  TOPSIS  method  we  propose.  For 
example,  assume  that  we  have  sensors  that  are  weather  dependent.  That  is,  as  a  storm  front  moves  into  range,  the 
strength  of  the  information  reported  by  various  sensors/nodes  can  deteriorate,  and  then  return  to  normal  after  the 
storm  passes.  This  can  be  modeled  by  using  weights  that  are  time  varying  in  NoRT. 

V.  TOPSIS  &  NoRT 

TOPSIS  was  first  discussed  in  ifTH.  Our  two  main  sources  are  given  above,  and  Il27l. 


1)  We  start  with  amxn  Decision  Matrix  D  =  dm,n,  where  the  rows  z  =  1, m  are  the  alternatives  (e.g.  nodes), 
and  the  columns  j  =  1,  ...,n  are  the  criteria  (e.g.  centrality  measures,  node  asset  values).  The  alternatives  are 
{.Ai},  and  the  criteria  are  {/Cj}.  The  set  of  criteria  that  is  indexed  hy  j,  consists  of  both  “henefit”  criteria 
and  “cost”  criteria  /C“.  If  there  are  no  cost  criteria,  the  TOPSIS  algorithm  simplifies.  For  NoRT  we  are  able 
to  restrict  to  only  benefit  criteria.  This  is  because  any  cost  effects  can  be  incorporated  into  the  non-centrality 
criteria  called  asset  criticality  (AS).  The  higher  the  AS  value  is  of  a  node  the  more  important  that  node  is  to 
us.  Therefore,  with  respect  to  NoRT,  {/Cj}  =  /C+,  and  JC~  =  0. 

Hence,  we  simplify  notation  and  use  A  to  denote  the  set  alternatives,  and  /C  to  denote  the  set  of  criteria. 

2)  We  next  form  a  normalized  decision  matrix  D  =  Sij  by  normalizing  every  entry  in  Dj,  the  jth  column  of 
D,  by  its  column  norm  \\Dj\\. 

8i,j  =  z  =  1,  =  1,  ...,n.  (6) 

/  m 


We  note  that  if  Dk  is  already  normalized  (as  in  the  Gould  vector  for  Eigenvalue  centrality)  that  5i^k  =  di^k- 
3)  We  now  assign  a  weight  of  0  <  wj  to  each  column  and  form  the  weighted  normalized  decision  matrix 

V  =  Vij 


Vij=wj-Sij  i  =  1,  =  1,  (7) 

Even  though  it  is  not  necessary  we  normalize  the  weights  so  that 
4)  We  next  calculate  the  positive  ideal  solution  A+  and  the  negative  ideal  solution  A~ . 

and  (8) 

where  (9) 

vf  =  maxvij  ,  and  (10) 

i 

v~  =  minvij  ,  j  =  l,...,n.  (11) 

i 


Thus,  A+  corresponds  to  finding  fhe  largest  value  in  every  criteria  (column),  and  A~  to  the  minimal. 

5)  Now  we  go  through  attribute  by  attribute  (row)  to  see  how  the  actual  values  compare  to  the  positive  ideal, 
and  the  negative  ideal.  We  calculate  the  Positive  (Negative)  Separation  (S~  )  between  each  alternative 
and  the  positive  ideal  solution  (negative  ideal  solution). 


n 


(12) 

(13) 


6) 


We  next  calculate  the  relative  closeness  of  alternate  z  to  the  ideal  positive  solution.  We  want  to  be  close  to 
the  ideal  positive  solution  and  far  from  the  ideal  negative  solution  (see  ll30l).  The  relative  closeness  of  the  z 
alternate  (node  for  us)  is  defined  as 


Q 


s- 

s-  +  s^  • 


(14) 


When  Ci  =  1  {Ci  =  0)  alfemafive  z  is  fhe  besf  (worsf)  solution,  that  is  it  coincides  with  A+  {A~). 

Note  that  if  we  have  a  weighting  wj  and  replace  it  with  K  ■  Wj  the  terms  and  S~  are  multiplied  by 
K,  but  the  Ci  values  are  unchanged.  In  particular  if  every  column  is  equi-weighted  at  1,  and  we  change 
the  normalization  to  Wj  =  1/n  the  Ci  are  unchanged.  We  choose  the  default  normalization  of  Wj  =  1  if  all 
weights  are  equal  so  that  we  can  ignore  step  (2)  above. 

7)  Now  we  rank  the  alternatives,  from  highest  to  lowest,  via  the  Ci  value. 
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V-A.  Example  2 

For  completeness,  we  include  the  following  example  for  algorithm  verification  against  a  published  example  of 
TOPSIS.  We  note  that  this  example  has  nothing  to  do  with  graphs  or  nodes.  The  alternatives  are  years  2008,  ..., 
2012;  i  =  1,...,5))  and  criteria  (various  rates  SR,PR,  and  CR,  j  =  1,2,3.  We  use  a  weighting  of  wj  =  1/3  for 
each  criteria.  We  see  that  the  rankings  of  the  years  in  descending  order  are  2010,  2012,  2011,  2008,  2009,  as  shown 
hy  our  excel  calculations  helow  (and  agreeing  with  0). 


YEAR  SR 

PR 

CR 

SRn=vil 

PRn=vi2 

CRn=vi3 

A+  =vj+  A-  =  vj- 

2008  0.950 

0.953 

0.950 

0.144 

0.148 

0.153 

0,152  0.144 

2009  1.000 

0.900 

0.902 

0.152 

0.140 

0.145 

0,153  0.140 

2010  0.974 

0.975 

0.946 

0.148 

0.152 

0.152 

0,153  0.145 

2011  0.984 

0.982 

0.903 

0.149 

0.153 

0.145 

2012  1.000 

0.974 

0.925 

0.152 

0.152 

0.149 

S+ 

S- 

Ci 

RANK 

S1+ 

0.009 

SI- 

0.011 

Cl 

0.561 

4,000 

S2+ 

0.015 

S2- 

0.008 

C2 

0.337 

5,000 

S3+ 

0.004 

S3- 

0.014 

C3 

0.773 

1,000 

S4+ 

0.008 

S4- 

0.014 

C4 

0.634 

3,000 

S5+ 

0.004 

S5- 

0.014 

C5 

0.772 

2,000 

Fig.  5:  Example  2  (Table  2)  from  0 
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V-B.  Return  to  Interesting  Example  1 

Once  again,  we  consider  the  Kite  as  shown  in  Fig.  This  time  we  will  use  the  normalizations  that  we  discussed 
at  the  beginning  of  the  paper,  and  we  will  add  a  column  for  asset  criticality  (AC).  We  start  with  a  default  value  of 
1  for  every  row  in  AC.  Of  course  this  node  ranking  agrees  with  the  ranking  done  without  any  asset  criticality. 


X 


• 

DC 

Ck>ser>e&s 

Between 

EC 

Squares  and  sums 

8 

S' 

C 

Beverty 

1 

0.4444 

0.059 

0.023 

1.00000 

0.197531 

0.X348 

0.00054 

1.00000 

0.3266 

0.3120 

0.0410 

0.3522 

0.68287 

0  44212 

0.39300 

Beverty 

6 

Andre 

2 

0.4444 

0.059 

0.023 

1.00000 

0.197531 

0.00348 

O.OOOS4 

1.00000 

0.3266 

0.3120 

0.0410 

0.3522 

0.68287 

0  44212 

0.39300 

Artdre 

6 

Cerot 

3 

0.3333 

0056 

0.000 

0.11155 

0.111111 

0.00314 

000000 

0.65861 

0.2449 

0  2962 

0.0000 

0.2858 

0.76035 

0  34007 

0.30904 

Carol 

8 

rernendo 

4 

0.SSS6 

0.071 

0.231 

1 12913 

0.308642 

0.00504 

0  05358 

1.27493 

0.4082 

0.3755 

0.4097 

0.3977 

0.30204 

067986 

0.69240 

Fernando 

1 

Garth 

5 

0.SSS6 

0.071 

0.231 

1.12913 

0.308642 

0.00504 

0.05358 

1.27493 

0.4082 

0.3755 

0.4097 

0.3977 

0.30204 

0.67986 

0.69240 

Garth 

1 

Ed 

6 

0.3333 

0  056 

0.000 

0.81155 

0.111111 

0.00314 

0  00000 

0.65861 

0.2449 

02962 

0.0000 

0.2858 

0.76035 

034007 

0.30904 

Ed 

8 

Diane 

7 

0.6667 

0  067 

0.102 

1.36572 

0444444 

0.00449 

001037 

1.86519 

0.4899 

0  3543 

0.1803 

o.aio 

0.50847 

0.67111 

0.56894 

Diane 

4 

Heather 

8 

0.3333 

0.067 

0.389 

0.55609 

0.111111 

0.00449 

0.15123 

0.30924 

0.2449 

0.3543 

0.6883 

0.1959 

0.37651 

0.75166 

0.66626 

Heather 

3 

Ike 

9 

02222 

0048 

0.222 

0.13649 

0.049383 

0.00230 

0  04938 

0.01863 

0.1633 

02539 

03933 

0.0481 

0.62922 

0  41013 

0.39460 

Ike 

5 

iane 

10 

o.ini 

0034 

0.000 

0.03169 

0.012346 

0.X116 

0.00000 

0.00100 

0.0816 

0.1798 

0.0000 

0.0112 

0.94840 

0.00000 

0.00000 

lane 

lO 

1.36083 

0.18909 

056500 

2.83922 

♦  Ideal 

0.4899 

0.3755 

0.6883 

0.4810 

Ideal 

0.0816 

0.1798 

0.0000 

0.0112 

Fig.  6:  Node  ranking  based  on  4  closeness  criteria 


A 

# 

DC 

Closeness 

Between 

EC 

AC 

Squares  and  sums 

R 

S+ 

S- 

C 

Beverly 

1 

0.4444 

0.059 

0.023 

1.00000 

1.00000 

0.197531 

0.00348 

0.00054 

1.00000 

1.00000 

0.3266 

0.3120  0,0410 

0.3522 

0.31622777 

0.68287 

0.44212 

0.39300 

Beverly 

L  ^ 

Andre 

2 

0.4444 

0.059 

0.023 

1.00000 

1.00000 

0.197531 

0.00348 

0.00054 

1.00000 

1.00000 

0.3266 

0.3120  0,0410 

0.3522 

0.31622777 

0.68287 

0.44212 

0.39300 

Andre 

6 

Carol 

3 

0.3333 

0.056 

0.000 

0.81155 

1.00000 

0.111111 

0.00314 

0.00000 

0,65861 

1.00000 

0.2449 

0.2962  0,0000 

0.2858 

0.31622777 

0.76035 

0.34007 

0.30904 

Carol 

8 

Fernando 

4 

0.5556 

0.071 

0.231 

1.12913 

1.00000 

0.308642 

0.00504 

0.05358 

1.27493 

l.OOOOOj 

0.4082 

0.3755  0,4097 

0.3977 

0.31622777 

0.30204 

0.67986 

0.69240 

Fernando 

1 

Garth 

5 

0.5556 

0.071 

0.231 

1.12913 

1.00000 

0.308642 

0.00504 

0.05358 

1,27493 

i.ooooo' 

0.4082 

0.3755  0,4097 

0.3977 

0.31622777 

0,30204 

0.67986 

0.69240 

Garth 

1 

Ed 

6 

0.3333 

0.056 

0.000 

0.81155 

1.00000 

0.111111 

0.00314 

0.00000 

0.65861 

1.00000 

0.2449 

0.2962  0,0000 

0.2858 

0.31622777 

0.76035 

0.34007 

0.30904 

Ed 

8 

Diane 

7 

0.6667 

0.067 

0.102 

1.36572 

1.00000 

0,444444 

0.00449 

0.01037 

1.86519 

1.00000 

0.4899 

0.3543  0,1803 

0.4810 

0.31622777 

0.50847 

0.67111 

0.56894 

Diane 

4 

Heather 

8 

0.3333 

0.067 

0.389 

0.55609 

1.00000 

0.111111 

0.00449 

0.15123 

0.30924 

1.00000 

0.2449 

0.3543  0.6883 

0.1959 

0.31622777 

0.37651 

0.75166 

0.66626 

Heather 

3 

Ike 

9 

0.2222 

0.048 

0.222 

0.13649 

1.00000 

0.049383 

0.00230 

0.04938 

0.01863 

I.ooooo] 

0.1633 

0.2539  0,3933 

0.0481 

0.31622777 

0.62922 

0.41013 

0.39460 

Ike 

5 

Jane 

10 

0.1111 

0.034 

0.000 

0.03169 

1.00000 

0.012346 

0.00116 

0.00000 

0.00100 

I.ooooo' 

0.0816 

0.1798  0.0000 

0.0112 

0.31622777 

0.94840 

0.00000 

0.00000 

Jane 

lO 

1.36083 

0.18909 

0.56500 

2.83922 

3.16228 

+  Ideal 

0.4899 

0.3755  0,6883 

0.4810 

0.31623 

-  Ideal 

0.0816 

0.1798  0,0000 

0.0112 

0.3162 

Fig.  7:  Node  ranking  based  on  4  closeness  criteria  and  all  equal  Asset  Critieria-of  course  there  is  no  difference. 
This  will  only  happen  if  the  AC  values  are  not  all  equal. 
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Node 

TOPSIS  Rank 

4 

first 

5 

first 

8 

third 

7 

fourth 

9 

fifth 

1 

sixth 

2 

sixth 

3 

eighth 

6 

eighth 

10 

tenth 

TABLE  II:  Rank 


Using  Table  II,  we  see  the  rank  of  each  node  (alternate)  based  on  DC,  CC,  BC,  or  EC.  When  we  use  TOPSIS 
with  equal  weighting  for  each  criteria  we  arrive  at  the  above  ranking,  which  does  not  agree  with  any  of  the  other 
four! 

Hence,  TOPSIS  can  be  used  successfully  when  we  have  to  balance  different  measures  of  centrality  and  node 
importance  against  each  other. 
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VI.  User  Manual 

This  program  uses  SageMath  and  Python  to  calculate  the  TOPSIS  of  a  graph  based  on  several  measures.  These 
measures  are:  Closeness  Centrality  (CC);  Betweenness  Centrality  (BC);  Degree  Centrality  (DC);  and  Eigenvector 
Centrality  (EC).  It  may  also  take  in  optional  Asset  Critically  (AC)  measurements.  AC  takes  in  a  separate  weights 
file  than  do  the  centrality  measurements.  Eor  the  AC  weights  file,  use  a  plain  editor  and  save  as  a  CSV  file.  The 
first  number  represents  the  node  number  and  should  be  an  integer.  There  should  be  a  comma  to  separate  the  first 
number  and  the  second  number.  The  second  number  represents  the  weight  and  should  be  a  floating  point  decimal. 
Eor  example  1  would  be  1.0  and  .1  would  be  0.1. 

The  script  takes  in  an  input  file  of  JSON  for  the  specifications  of  the  graph.  Alternatively,  it  also  accepts  a  TSV 
file  defining  the  edges  of  the  graph.  In  both  cases,  nodes  are  defined  as  integers  greater  than  zero.  In  the  following 
steps,  replace  any  words  inside  the  brackets  ([.  .  .])  with  your  (the  user’s)  information.  This  program  is  meant  to 
be  run  on  a  linux  machine. 

1)  To  install  in  SageMath  Python  environment,  go  to  command  line  and  type:  sage  -python  -m  pip 
install  pathlib  argparse 

2)  To  print  a  short  description  of  all  command  line  options  for  the  script,  type:  sage  -python  [program 
name] .py  -h 

3)  The  standard  line  to  run  the  script  using  a  JSON  file  is  as  follows: 

sage  -python  [program  name] .py  [file  name] . js  -[optional  measure]  [optional 
weight ] 

4)  To  run  the  script  using  a  TSV  file,  replace  .js  with  .tsv 

5)  To  add  an  optional  argument  for  asset  criticality  (AC),  type:  sage  -python  [program  name]  .py 
[file  name] . js  -[optional  measure]  [optional  weight]  -AC  [file  name] . csv 

6)  The  output  of  your  program  should  create  a  text  file  which  contains  data  similar  to  this: 

Rankings  for  graph  kite.js 
Number  TOPSIS  Closeness  C  Value 
8  1.000000 

4  0.595238 

5  0.595238 

9  0.571429 

7  0.261905 

1  0.059524 

2  0.059524 

3  0.000000 

6  0.000000 

10  0.000000 

TAB  EE  3:  Rankings  for  BC 

sage  -python  topsisRK.py  kite.js  -BC  1.0 

VII.  Experiment 

Once  again  looking  at  the  Kite  graph  from  Eigure  3,  we  combine  the  four  different  centrality  measures  into 
combinations  of  pairs.  Eirst,  we  initialize  one  of  the  measures  to  0.1  while  keeping  the  other  at  1.0.  Then  we 
shuffle  which  measurements  are  paired  with  each  other  and  repeat  these  steps  for  each  pair. 

Below  are  a  series  of  tables  that  show  a  few  results  of  the  different  centrality  measurement  pairs  and  their 
respective  weights.  Each  table  contains  the  node  ranking,  node  number,  and  centrality  values. 


Rank 

1 

2 

2 

4 

5 

6 
6 
8 
8 
8 
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Rank  Node  Number  TOPSIS  Closeness  C  Value 


1 

8 

0.965644 

2 

4 

0.595898 

2 

5 

0.595898 

4 

9 

0.569989 

5 

7 

0.266768 

6 

1 

0.068656 

6 

2 

0.068656 

8 

3 

0.023161 

8 

6 

0.023161 

10 

10 

0.000000 

TABEE  4:  Rankings  of  BC 

=  1.0,  DC  = 

0.1 

-python  topsisRK.py 

kite . j  s  -BC 

1.0  -E 

Rank 

Node  Number  TOPSIS  Closeness  C  Value 

1 

7 

0.889427 

2 

4 

0.792331 

2 

5 

0.792331 

4 

1 

0.582399 

4 

2 

0.582399 

6 

8 

0.419774 

7 

3 

0.390916 

7 

6 

0.390916 

9 

9 

0.216527 

10 

10 

0.000000 

TABEE  5:  Rankings  of  BC 

=  0.1,  DC  = 

1.0 

-python  topsisRK.py 

kite . j  s  -BC 

0.1  -E 

Rank 

Node  Number  TOPSIS  Closeness  C  Value 

1 

8 

0.960232 

2 

4 

0.596197 

2 

5 

0.596197 

4 

9 

0.568828 

5 

7 

0.268307 

6 

1 

0.076071 

6 

2 

0.076071 

8 

3 

0.038360 

8 

6 

0.038360 

10 

10 

0.000000 

TABLE  6:  Rankings  of  BC  =  1.0,  EC  =  0.1 

sage  -python  topsisRK.py  kite.js  -BC  1.0  -EC  0.1 
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Rank 

Node  Number 

TOPSIS  Closeness  C  Value 

1 

7 

0.902491 

2 

4 

0.815627 

2 

5 

0.815627 

4 

1 

0.702901 

4 

2 

0.702901 

6 

3 

0.570287 

6 

6 

0.570287 

8 

8 

0.408707 

9 

9 

0.110554 

10 

10 

0.000000 

TABLE  7:  Rankings  of  BC  =  0.1,  EC  =  1.0 

sage  -python  topsisRK.py  kite.js  -BC  0.1  -EC  1.0 


According  to  Tables  4  &  5,  we  observe  that  ranking  varies  depending  on  which  centrality  measures  we  use. 
When  BC  has  a  greater  weight  than  DC,  node  8  is  ranked  first.  On  the  other  hand,  if  DC  has  a  greater  weight, 
then  node  7  is  ranked  much  higher  than  node  8. 

Erom  these  results  we  see  that  the  ranking  of  a  node  is  not  only  affected  by  its  centrality  but  also  by  its  position 
in  the  network  structure.  Node  8  would  have  the  highest  BC  in  the  Kite  graph  because  it  is  placed  in  a  position 
connected  to  the  greatest  number  of  nodes  that  are  not  connected  to  each  other.  Conversely,  node  7  would  be  ranked 
much  higher  in  than  node  8  in  terms  of  DC  because  it  has  more  alternative  routes  than  other  nodes  to  reach  its 
goal. 

Tables  5  &  7  show  that  node  7  consistently  has  a  high  ranking  for  both  DC  and  EC.  As  stated  before,  DC  is  the 
number  of  links  connected  to  a  node  and  the  basic  concept  behind  EC  is  that  an  important  node  is  connected  to 
important  neighbors.  Given  the  position  of  node  7  in  the  Kite  graph,  it  would  make  sense  for  it  to  be  ranked  first 
for  both  centrality  measures.  Note  that  even  though  EC  and  DC  rank  the  nodes  the  same,  the  TOPSIS  Closeness 
C  values  for  them  are  normalized  differently  (think  of  100  being  higher  than  15,  and  20  being  higher  than  15). 

Eurther  testing  is  done  by  using  asset  criticality  (AC).  We  start  by  taking  the  different  centrality  measures  and 
finding  the  node  rank.  Then  we  take  the  highest  ranking  node  and  change  its  asset  criticality  to  0.1.  The  lowest 
ranking  node  is  given  a  higher  asset  criticality  of  3.0.  The  other  nodes  are  kept  constant  with  an  asset  criticality  of 
1.0.  Below  are  a  series  of  tables  that  show  the  rankings  of  each  centrality  measure  after  asset  criticality  has  been 
changed  between  high  and  low  ranked  nodes. 


Rank 

Node  Number 

AC  Weights  =  2nd  column  weightsRK.csv 

1 

10 

3.0 

2 

8 

0.1 

3 

4 

1.0 

3 

5 

1.0 

5 

9 

1.0 

6 

7 

1.0 

7 

1 

1.0 

7 

2 

1.0 

9 

3 

1.0 

9 

6 

1.0 

TAB  EE  8:  Switched  AC  weights  for  BC 

sage  -python  topsisRK.py  kite.js  -BC  1.0  -AC  weightsRK . csv 


15 


Rank  Node  Number  AC  Weights  =  2nd  column  weightsRK.csv 


1 

10 

3.0 

2 

4 

1.0 

2 

5 

1.0 

4 

1 

1.0 

4 

2 

1.0 

6 

7 

0.1 

7 

3 

1.0 

7 

6 

1.0 

7 

8 

1.0 

10 

9 

1.0 

TAB  EE  9:  Switched  AC  weights  for  DC 

-python 

topsisRK . py 

kite.js  -DC  1.0  -AC  weightsR 

Rank 

Node  Number  AC  Weights  =  2nd  column  weightsRK.csv 

1 

10 

3.0 

2 

5 

1.0 

3 

7 

1.0 

3 

8 

1.0 

5 

1 

1.0 

5 

2 

1.0 

7 

3 

1.0 

7 

6 

1.0 

9 

9 

1.0 

10 

4 

0.1 

TABEE  10:  Switched  AC  weights  for  CC 

-python 

topsisRK . py 

kite.js  -CC  1.0  -AC  weightsR 

Rank 

Node  Number  AC  Weights  =  2nd  column  weightsRK.csv 

1 

10 

3.0 

2 

4 

1.0 

2 

5 

1.0 

4 

1 

1.0 

4 

2 

1.0 

6 

3 

1.0 

6 

6 

1.0 

8 

7 

0.1 

9 

8 

1.0 

10 

9 

1.0 

TABLE  1 1 :  Switched  AC  weights  for  EC 

sage  -python  topsisRK.py  kite.js  -EC  1.0  -AC  weightsRK.csv 

Eooking  at  Table  11,  we  see  the  nodes  behave  as  expected.  When  we  change  asset  criticality,  the  lowest  ranking 
node  for  each  centrality  measure  becomes  the  highest  ranking  node.  The  nodes  which  were  ranked  between  the 
highest  and  lowest  nodes  stayed  in  the  middle,  and  the  previous  first  ranked  node  gets  ranked  much  lower. 


VII-A.  Sagemath  Python  Code:  topsisRK.py 
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0)  C  01  O 

0)  01  01  01 

iM  CO  3  > 

Cn  o  4J  C 

01  rH  01  01 
T5  O  tJi 
I  I  I 
>.  >1  >1  01 

4->  4J  4J  I 


cn  a  T5  T3 


tj'  x: 
0!  O- 
CO  0! 


01  x:  CO  o 
I  ^  H  c: 
t  CO  vj 


fill 

lO  cO  cO  cO 


^3  01  C  01  > 


CO  O  4J  O  lO 


01  Co  O  4-)  CTI 


TO  T3  TO  T3  t#’ 


01  01  01  01  01 


CO  CO  c:  01  G 


CO  O  01 
CO  x:  > 

CO  G)  4-)  CO  • 


Q  O  O  4-)  CO  70 


M  O  a  -H 
o  a  G 
as  CO 


CO  G  M  4-1 

G  05  X2  -H  >1  CO 

-G  Q-  th  G  a  • 

G  CO  M  •  S  >1 

4J  G  01  G  a 

CO  G  4-1  CO  G  -H 
G  G  U 


OOOOOOSSOS 


O  O  O  O  O 
U  U  U  U  O 


01  01  01  01  01 


DO 

01 

Co  w 
G  DC 
Cd  05  4-> 
DC  DC  G 
Cd  Cd  4J  4J  01 

DC  DC  05  05  6 

>1  4J  O  O  S 

01  G  M  M  O 


05  G 
01  01 
s  s 

01  G 


■»-CMCO-^ir)CDh'C0030'«-C 


#  computes  a  centrality  measurement  of  the  nodes  in  graph  (g)  with  eigenvectors 
def  centrality_eigenvector (g) : 

print  "  calculating  eigenvectors..." 


fK  =  float (k) 

elif  isinst ance (k, AlgebraicNumber ) : 
sK  =  str  (k) 

if  floatRangeRE .match (sK) : 
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>  O  01  ‘ 
ca  o  a  t 

H  -H  C  fO 
(J)  Q  -H  (JJ 
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e  o  -o 

IT5  g  c 

S  fO  fO 

o;  z 


)  cn  T5 

i  C  O 
:  -H  c 


dJ  o  j<j  0- 


(U  U  -H  u 

■— I  -H  [14  -H 

•H  Q  4J  Q 
■P  a  04  04 


>  -H  CC  . 

cn  4J  X  o  T 

u  (T5  Q)  Cd 


a)0404^>i.i-’M>i04M 


G  dJ  4J  0 


d)  ^  G  G  • 


T^locDr^coo)0^-<Nco•^u^cD^^coo>o■^-(NcoT^lncDr^ooCT^Ol-<Nco■^lncD^^coo50■^-cvJco■>J■lncDl^coo>o■^-c\J(r)■^lncDr^coc350^-OJco•^ln<D^>-coa>o■^-c\JcoT^l^) 

■^■^■^■>j-->j-->j-inininu^ininu^u^inu^cDcDcDcDcDcDcDcDcD<Dr^h'h'h'h'h~h'h'r^r^cococooooococoooooooo>05CT>cT>a>05C35C35C35C350000000000-»-i--»-^^^ 
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C  0-0- 

iB  -H  C  C 

e  e 

iB  ni  II  II 


M  q  q  q 

Dj  -h  -h  -h  q 

M  M  M  O 

*  0-  0-  0-  9-1 


cDr^coo)0^-cvico•^ln<Dl^cocT>o■^-(NcOT^LncDr^coo)Ol-<Nco•^ln(D^'C0050■^-(Nco■>J■lr)cD^-ooa^o■^-<^Jco■^l^)cD^-coc350■^-<Nco■>J■lr)<D^>-coo>o 

■«-'«-'«-T-cvjojcvjC'j<N<N<N<NC\jc\jco(r)cocococococococo->j--^-^-^-^-^'^'^'^'^Lnioinioinininif5inin<D<D<DcDcDcDcDcDcDcDr^r^r^r^r^r^i^i^i^i^co 

(N(N(NCViCVJCViCviCvi<NC\l<NC\IC\JC\JC\JC\J(N(N(NC\J(NCViCVJCVJCVJ<NC\l<N<N<NC\JC\JC\JC\JC\J(N(N(NCVJCViCVJCVJCVJC\l<NC\l<NC\JC\JC\JC\J(N(N(N(N(NCViCviCviCviC\l<NC\l<NW 
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